Conformational and dynamic properties of the KH1 domain of FMRP and its fragile X syndrome linked G266E variant.

The Fragile X messenger ribonucleoprotein (FMRP) is a complex, multi-domain protein involved in interactions with various macromolecules, including proteins and coding/non-coding RNAs. The three KH domains (KH0, KH1 and KH2) within FMRP are recognized for their roles in mRNA binding. In the context of Fragile X syndrome (FXS), over-and-above CGG triplet repeats expansion, three specific point mutations have been identified, each affecting one of the three KH domains (R138QKH0, G266EKH1, and I304NKH2) resulting in the expression of non-functional FMRP. This study aims to elucidate the molecular mechanism underlying the loss of function associated with the G266EKH1 pathological variant. We investigate the conformational and dynamical properties of the isolated KH1 domain and the two KH1 site-directed mutants G266EKH1 and G266AKH1. Employing a combined in vitro and in silico approach, we reveal that the G266EKH1 variant lacks the characteristic features of a folded domain. This observation provides an explanation for functional impairment observed in FMRP carrying the G266E mutation within the KH1 domain, as it renders the domain unable to fold properly. Molecular Dynamics simulations suggest a pivotal role for residue 266 in regulating the structural stability of the KH domains, primarily through stabilizing the α-helices of the domain. Overall, these findings enhance our comprehension of the molecular basis for the dysfunction associated with the G266EKH1 variant in FMRP.

Unveiling the Cleavage Mechanism of an RNA Model Compound on the whole pH Scale: Computations Meet Experiments in the Determination of Reaction Rates.

The knowledge of the mechanism of reactions occurring in solution is a primary research line both in the context of theoretical-computational chemistry and in the field of organic and bio-organic chemistry. Given the importance of the hydrolysis of nucleic acids in life-related phenomena, here we present a combined experimental and computational study on the cleavage of an RNA model compound. This phosphodiester features a cleavage rate strictly dependent on the pH with three different dependence domains. Such experimental evidence, highlighted by an in-depth kinetic investigation, unequivocally suggests a change in the reaction mechanism along the pH scale. In order to interpret the data and to explain the experimental behavior, we have applied a theoretical-computational procedure, involving a hybrid quantum/classical approach, able to model chemical reactions in complex environments, i.e. in solution. This study turns out to quantitatively reproduce the experimental data with accuracy and, in addition, provides useful mechanistic insight into the transesterification process of the investigated compound. The study indicates that the cleavage can occur through an ANDN, an AN+DN, and a DNAN mechanism depending on the pH values.

Modeling the temperature dependence of the fluorescence properties of Indole in aqueous solution.

In a recent paper, we proposed a scheme to describe the relaxation mechanism of the excited Indole in aqueous solution, involving the fluctuations among the diabatic electronic states 1Lb, 1La and 1πσ∗. Such a theoretical and computational model reproduced accurately the available experimental data at room temperature. Following these results, in the present work, we model the complex temperature dependence of the fluorescence properties of Indole in aqueous solution, with results further validating the proposed relaxation scheme. This scheme is able to explain the temperature effects on the fluorescence behavior indicating the water fluctuations as the main cause of (i) the stabilization of the dark state (1πσ∗) and (ii) the increase in temperature of the kinetics of the irreversible transition towards such a state.

A Theoretical-Computational Study of Phosphodiester Bond Cleavage Kinetics as a Function of the Temperature.

The hydrolysis of the phosphodiester bond is an important chemical reaction involved in several biological processes. Here, we study the cleavage of this bond by means of a theoretical-computational method in a model system, the dineopentyl phosphate. By such an approach, we reconstructed the kinetics and related thermodynamics of this chemical reaction along an isochore. In particular, we evaluated the kinetic constants of all the reaction steps within a wide range of temperatures, mostly corresponding to conditions where no experimental measures are available due to the extremely slow kinetics. Our results, in good agreement with the experimental data, show the robustness of our theoretical-computational methodology which can be easily extended to more complex systems.

A partial human LCK defect causes a T cell immunodeficiency with intestinal inflammation.

Lymphocyte-specific protein tyrosine kinase (LCK) is essential for T cell antigen receptor (TCR)-mediated signal transduction. Here, we report two siblings homozygous for a novel LCK variant (c.1318C>T; P440S) characterized by T cell lymphopenia with skewed memory phenotype, infant-onset recurrent infections, failure to thrive, and protracted diarrhea. The patients' T cells show residual TCR signal transduction and proliferation following anti-CD3/CD28 and phytohemagglutinin (PHA) stimulation. We demonstrate in mouse models that complete (Lck-/-) versus partial (LckP440S/P440S) loss-of-function LCK causes disease with differing phenotypes. While both Lck-/- and LckP440S/P440S mice exhibit arrested thymic T cell development and profound T cell lymphopenia, only LckP440S/P440S mice show residual T cell proliferation, cytokine production, and intestinal inflammation. Furthermore, the intestinal disease in the LckP440S/P440S mice is prevented by CD4+ T cell depletion or regulatory T cell transfer. These findings demonstrate that P440S LCK spares sufficient T cell function to allow the maturation of some conventional T cells but not regulatory T cells-leading to intestinal inflammation.

Molecular Modeling of the Deamidation Reaction in Solution: A Theoretical-Computational Study.

In this work, a theoretical-computational method is applied to study the deamidation reaction, a critical post-translational modification in proteins, using a simple model molecule in solution. The method allows one to comprehensively address the environmental effect, thereby enabling one to accurately derive the kinetic rate constants for the three main steps of the deamidation process. The results presented, in rather good agreement with the available experimental data, underline the necessity for a rigorous treatment of environmental factors and a precise kinetic model to correctly assess the overall kinetics of the deamidation reaction.

ERAP1 and ERAP2 Haplotypes Influence Suboptimal HLA-B*27:05-Restricted Anti-Viral CD8+ T Cell Responses Cross-Reactive to Self-Epitopes.

The human leukocyte antigen (HLA)-B*27 family of alleles is strongly associated with ankylosing spondylitis (AS), a chronic inflammatory disorder affecting the axial and peripheral joints, yet some HLA-B*27 variants not associated with AS have been shown. Since no major differences in the ligandome of associated compared to not-associated alleles have emerged, a plausible hypothesis is that the quantity rather than the quality of the presented epitopes makes the difference. In addition, the Endoplasmic Reticulum AminoPeptidases (ERAPs) 1 and 2, playing a crucial role in shaping the HLA class I epitopes, act as strong AS susceptibility factors, suggesting that an altered peptidome might be responsible for the activation of pathogenic CD8+ T cells. In this context, we have previously singled out a B*27:05-restricted CD8+ T cell response against pEBNA3A (RPPIFIRRL), an EBV peptide lacking the B*27 classic binding motif. Here, we show that a specific ERAP1/2 haplotype negatively correlates with such response in B*27:05 subjects. Moreover, we prove that the B*27:05 allele successfully presents peptides with the same suboptimal N-terminal RP motif, including the self-peptide, pDYNEIN (RPPIFGDFL). Overall, this study underscores the cooperation between the HLA-B*27 and ERAP1/2 allelic variants in defining CD8+ T cell reactivity to suboptimal viral and self-B*27 peptides and prompts further investigation of the B*27:05 peptidome composition.

Unveiling the Excited State Dynamics of Indole in Solution.

In this paper, we reconstruct in detail the dynamics of the emitting electronic excited state of aqueous indole, investigating its relaxation mechanism and kinetics to be related to the time-dependent fluorescence signal. Taking advantage of the results shown in a very recent paper, we were able to model the relaxation process in solution in terms of the transitions between two gas-phase singlet electronic states (1La and 1Lb), subsequently irreversibly relaxing to the gas-phase singlet dark state (1πσ*). A comparison of the results with the available experimental data shows that the relaxation mechanism we obtain by our theoretical-computational model is reliable, reproducing rather accurately all the experimental observables.

A study of cyanidin/alginate complexation: Influence of pH in assembly and chiral properties.

Cyanidin 3-O-glucoside (CND) is a frequently-used anthocyanin that has excellent antioxidant properties but a limited bioavailability in bloodstream. Complexation of CND with alginate can improve its therapeutic outcome. Here we have studied the complexation of CND with alginate under a range of pH values from 2.5 to 5. CND is positively charged at low pH, and becomes neutral, and then negatively charged as pH increases. CND/alginate complexation was studied by dynamic light scattering, transmission electron microscopy, small angle X-ray scattering, STEM, UV-Vis spectroscopy and circular dichroism (CD). CND/alginate complexes at pH 4.0 and 5.0 form chiral fibres with a fractal structure. At these pH values, CD spectra show very intense bands, which are inverted compared with free CND. Complexation at lower pH results in disordered polymer structures and CD spectra show the same features as for CND in solution. Molecular dynamics simulations suggest the formation of parallel CND dimers through complexation with alginate at pH 3.0, while at pH 4.0 CND dimers form in a cross like arrangement.

Modelling Complex Bimolecular Reactions in a Condensed Phase: The Case of Phosphodiester Hydrolysis.

(1) Background: the theoretical modelling of reactions occurring in liquid phase is a research line of primary importance both in theoretical-computational chemistry and in the context of organic and biological chemistry. Here we present the modelling of the kinetics of the hydroxide-promoted hydrolysis of phosphoric diesters. (2) Method: the theoretical-computational procedure involves a hybrid quantum/classical approach based on the perturbed matrix method (PMM) in conjunction with molecular mechanics. (3) Results: the presented study reproduces the experimental data both in the rate constants and in the mechanistic aspects (C-O bond vs. O-P bond reactivity). The study suggests that the basic hydrolysis of phosphodiesters occurs through a concerted ANDN mechanism, with no formation of penta-coordinated species as reaction intermediates. (4) Conclusions: the presented approach, despite the approximations, is potentially applicable to a large number of bimolecular transformations in solution and therefore leads the way to a fast and general method to predict the rate constants and reactivities/selectivities in complex environments.

Effect of Salts on the Conformational Dynamics of the Cytochrome P450 OleP.

Cytochrome P450 OleP catalytic activity is strongly influenced by its structural dynamic conformational behavior. Here, we combine equilibrium-binding experiments with all-atom molecular dynamics simulations to clarify how different environments affect OleP conformational equilibrium between the open and the closed-catalytic competent-forms. Our data clearly show that at high-ionic strength conditions, the closed form is favored, and, very interestingly, different mechanisms, depending on the chemistry of the cations, can be used to rationalize such an effect.

PyMM: An Open-Source Python Program for QM/MM Simulations Based on the Perturbed Matrix Method.

Quantum mechanical/molecular mechanics (QM/MM) methods are important tools in molecular modeling as they are able to couple an extended phase space sampling with an accurate description of the electronic properties of the system. Here, we describe a Python software package, called PyMM, which has been developed to apply a QM/MM approach, the perturbed matrix method, in a simple and efficient way. PyMM requires a classical atomic trajectory of the whole system and a set of unperturbed electronic properties of the ground and electronic excited states. The software output includes a set of the most common perturbed properties, such as the electronic excitation energies and the transitions dipole moments, as well as the eigenvectors describing the perturbed electronic states, which can be then used to estimate whatever electronic property. The software is composed of a simple and complete command-line interface, a set of internal input validation, and three main analyses focusing on (i) the perturbed eigenvector behavior, (ii) the calculation of the electronic absorption spectrum, and (iii) the estimation of the free energy differences along a reaction coordinate.

A Simplified Treatment for Efficiently Modeling the Spectral Signal of Vibronic Transitions: Application to Aqueous Indole.

In this paper, we introduce specific approximations to simplify the vibronic treatment in modeling absorption and emission spectra, allowing us to include a huge number of vibronic transitions in the calculations. Implementation of such a simplified vibronic treatment within our general approach for modelling vibronic spectra, based on molecular dynamics simulations and the perturbed matrix method, provided a quantitative reproduction of the absorption and emission spectra of aqueous indole with higher accuracy than the one obtained when using the existing vibronic treatment. Such results, showing the reliability of the approximations employed, indicate that the proposed method can be a very efficient and accurate tool for computational spectroscopy.

Theoretical Evaluation of Sulfur-Based Reactions as a Model for Biological Antioxidant Defense.

Sulfur-containing amino acids, Methionine (Met) and Cysteine (Cys), are very susceptible to Reactive Oxygen Species (ROS). Therefore, sulfur-based reactions regulate many biological processes, playing a key role in maintaining cellular redox homeostasis and modulating intracellular signaling cascades. In oxidative conditions, Met acts as a ROS scavenger, through Met sulfoxide formation, while thiol/disulfide interchange reactions take place between Cys residues as a response to many environmental stimuli. In this work, we apply a QM/MM theoretical-computational approach, which combines quantum-mechanical calculations with classical molecular dynamics simulations to estimate the free energy profile for the above-mentioned reactions in solution. The results obtained, in good agreement with experimental data, show the validity of our approach in modeling sulfur-based reactions, enabling us to study these mechanisms in more complex biological systems.

Modeling Charge Transfer Reactions by Hopping between Electronic Ground State Minima: Application to Hole Transfer between DNA Bases.

In this paper, we extend the previously described general model for charge transfer reactions, introducing specific changes to treat the hopping between energy minima of the electronic ground state (i.e., transitions between the corresponding vibrational ground states). We applied the theoretical-computational model to the charge transfer reactions in DNA molecules which still represent a challenge for a rational full understanding of their mechanism. Results show that the presented model can provide a valid, relatively simple, approach to quantitatively study such reactions shedding light on several important aspects of the reaction mechanism.

Role of the membrane anchor in the regulation of Lck activity.

Theoretical work suggests that collective spatiotemporal behavior of integral membrane proteins should be modulated by boundary lipids sheathing their membrane anchors. Here, we show evidence for this prediction while investigating the mechanism for maintaining a steady amount of the active form of integral membrane protein Lck kinase (LckA) by Lck trans-autophosphorylation regulated by the phosphatase CD45. We used super-resolution microscopy, flow cytometry, and pharmacological and genetic perturbation to gain insight into the spatiotemporal context of this process. We found that LckA is generated exclusively at the plasma membrane, where CD45 maintains it in a ceaseless dynamic equilibrium with its unphosphorylated precursor. Steady LckA shows linear dependence, after an initial threshold, over a considerable range of Lck expression levels. This behavior fits a phenomenological model of trans-autophosphorylation that becomes more efficient with increasing LckA. We then challenged steady LckA formation by genetically swapping the Lck membrane anchor with structurally divergent ones, such as that of Src or the transmembrane domains of LAT, CD4, palmitoylation-defective CD4 and CD45 that were expected to drastically modify Lck boundary lipids. We observed small but significant changes in LckA generation, except for the CD45 transmembrane domain that drastically reduced LckA due to its excessive lateral proximity to CD45. Comprehensively, LckA formation and maintenance can be best explained by lipid bilayer critical density fluctuations rather than liquid-ordered phase-separated nanodomains, as previously thought, with "like/unlike" boundary lipids driving dynamical proximity and remoteness of Lck with itself and with CD45.

Complexation and organization of doxorubicin on polystyrene sulfonate chains: impacts on doxorubicin dimerization and quenching.

Anthracycline doxorubicin hydrochloride (DX) is a positively charged fluorescent drug, which in water self-associates into non-fluorescent antiparallel dimers upon increasing concentration and/or ionic strength. The positive charge of DX allows for complexation with negatively charged polymers and drug carriers. The fluorescence of DX following complexation with polyanion polystyrene sulfonate (PSS) is studied here. The fluorescence emission of DX decreases in the presence of PSS, being almost completely quenched when the ratio (R) of PSS monomers-to-DX molecules is larger than 10. Increasing R values over 30 results in a progressive recovery of fluorescence. The circular dichroism of PSS-DX complexes shows inverted characteristic bands of DX dimers suggesting the presence of parallel dimers at a concentration of DX below dimerization in water. Molecular dynamics studies corroborate a preferential orientation of DX into parallel dimers when interacting with PSS and show that DX molecules interact with a binding pocket of PSS monomers rather than with one single monomer. Increasing the ionic strength results in a recovery of fluorescence without an apparent release of DX from the PSS-DX complex as shown by DOSY NMR. PSS acts as a template for concentrating DX, triggering dimerisation and orienting DX molecules with their charged groups facing the negatively charged PSS monomers.

Effects of Environmental and Electric Perturbations on the pKa of Thioredoxin Cysteine 35: A Computational Study.

Here we present a theoretical-computational study dealing with the evaluation of the pKa of the Cysteine residues in Thioredoxin (TRX) and in its complex with the Thioredoxin-interacting protein (TXNIP). The free energy differences between the anionic and neutral form of the Cysteine 32 and 35 have been evaluated by means of the Perturbed Matrix Method with classical perturbations due to both the environment and an exogenous electric field as provided by Molecular Dynamics (MD) simulations. The evaluation of the free energies allowed us to show that the effect of the perturbing terms is to lower the pKa of Cysteine 32 and Cysteine 35 with respect to the free amino-acid. On the other hand, in the complex TRX-TXNIP, our data show an enhanced stabilization of the neutral reduced form of Cys 35. These results suggest that external electric stimuli higher than 0.02 V/nm can modulate the Cysteine pKa, which can be connected to the tight regulation of the TRX acting as an antioxidant agent.

Atomic-Level View of the Functional Transition in Vertebrate Hemoglobins: The Case of Antarctic Fish Hbs.

Tetrameric hemoglobins (Hbs) are prototypal systems for studies aimed at unveiling basic structure-function relationships as well as investigating the molecular/structural basis of adaptation of living organisms to extreme conditions. However, a chronological analysis of decade-long studies conducted on Hbs is illuminating on the difficulties associated with the attempts of gaining functional insights from static structures. Here, we applied molecular dynamics (MD) simulations to explore the functional transition from the T to the R state of the hemoglobin of the Antarctic fish Trematomus bernacchii (HbTb). Our study clearly demonstrates the ability of the MD technique to accurately describe the transition of HbTb from the T to R-like states, as shown by a number of global and local structural indicators. A comparative analysis of the structural states that HbTb assumes in the simulations with those detected in previous MD analyses conducted on HbA (human Hb) highlights interesting analogies (similarity of the transition pathway) and differences (distinct population of intermediate states). In particular, the ability of HbTb to significantly populate intermediate states along the functional pathway explains the observed propensity of this protein to assume these structures in the crystalline state. It also explains some functional data reported on the protein that indicate the occurrence of other functional states in addition to the canonical R and T ones. These findings are in line with the emerging idea that the classical two-state view underlying tetrameric Hb functionality is probably an oversimplification and that other structural states play important roles in these proteins. The ability of MD simulations to accurately describe the functional pathway in tetrameric Hbs suggests that this approach may be effectively applied to unravel the molecular and structural basis of Hbs exhibiting peculiar functional properties as a consequence of the environmental adaptation of the host organism.

Rationalizing Sequence and Conformational Effects on the Guanine Oxidation in Different DNA Conformations.

The effect of the environment on the guanine redox potential is studied by means of a theoretical-computational approach. Our data, in agreement with previous experimental findings, clearly show that the presence of consecutive guanine bases in both single- and double-stranded DNA oligomers lowers their reduction potential. Such an effect is even more marked when a G-rich quadruplex is considered, where the oxidized form of guanine is particularly stabilized. To the best of our knowledge, this is the first computational study reporting on a quantitative estimate of the dependence of the guanine redox potential on sequence and conformational effects in complex DNA molecules, ranging from single-stranded DNA to G-quadruplex.